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The energy content of large-scale jets in FR I radio galaxies is still an open issue. Here we show that upper 
limits on the high-energy and very high-energy 7-ray emission of the kpc-scale jet in M 87 radio galaxy im- 
posed by EGRET, Whipple, and — most importantly — HEGRA and HESS observations, provide important 
constraints on the magnetic field strength in this object. In particular, a non-detection of 7-ray radiation from 
the brightest part of this jet (knot A), expected from the inverse-Compton scattering of the starlight photons 
by the synchrotron-emitting jet electrons, implies that the magnetic field cannot be smaller than the cquipar- 
tition value (referring solely to the radiating ultrarelativistic electrons), and most likely, is even stronger. In 
this context, we point out several consequences of the obtained result for the large-scale jet structures in FR 
I radio galaxies and the M 87 jet in particular. For example, we discuss a potential need for amplification of 
the magnetic field energy flux along these jets (from sub-pc to kpc scales), suggesting a turbulent dynamo as a 
plausible process responsible for the aforementioned amplification. 
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1. INTRODUCTION 

Stawarz et al. (2003) considered very high en- 
ergy (VHE) 7-ray emission produced by the kpc-scale 
jets in nearby low-power radio galaxies of the FR I 
type [Tt| . Optical and X-ray emission detected re- 
cently from a number of such objects indicate that 
these jets are still relativistic on the kpc scale, and 
that they contain ultrarelativistic electrons with ener- 
gies up to 100 TeV (see a discussion in [13 and ref- 
erences therein). Therefore, some of the nearby FR I 
jets can be in principle VHE 7-ray emitters due to the 
inverse-Compton (IC) scattering off ambient photon 
fields which, at kpc distances from the active nuclei, 
are expected to be still relatively high. For example, 
following 19j it can be found that the bolometric en- 
ergy density of the stellar emission at 1 kpc from the 
center of a typical luminous elliptical galaxy is on av- 
erage f/star.bol ~ 10~ 9 er g cm " 3 . In the particular 
case of radio galaxy M 87, |17j show that comptoniza- 
tion of such a starlight radiation (that dominates over 
the energy densities of the other photon fields in the 
jet comoving frame, in particular over the energy den- 
sity of the synchrotron photons) within the kpc-scale 
jet (its brightest knot A) by the synchrotron-emitting 
electrons in the equipartition magnetic field can pos- 
sibly account for the TeV emission detected by the 
HEGRA Cherenkov Telescope from the direction of 
that source @ However, subse quen t observations of 
M 87 by the Whipple Telescope [H4 gave only upper 
limits for its emission in the 0.4 — 4 TeV photon en- 
ergy range, suggesting, although not strictly implying, 
variability of the VHE 7-ray signal. Such variability, 
clearly confirmed by the most recent HESS observa- 
tions — which established the presence of a variable 
source of VHE 7-ray emission within 0.1 deg (~ 30 
kpc) of the M 87 central region Q — questions the 



possibility that the kpc-scale jet is responsible for the 
3 - 4ct detections of M 87 by HEGRA and HESS tele- 
scopes. On the other hand, the upper limits imposed 
in this way put interesting constraints on the magnetic 
field within the M 87 large-scale jet, an issue which is 
of general importance in astrophysics. 

Here we comment on the high energy 7-ray emis- 
sion of the M 87 kpc-scale jet, resulting from the IC 
scattering on the stellar photon field. We take into 
account a relativistic bulk velocity of the emitting re- 
gion as well as Klein-Nishina regime of the electron- 
photon interaction. We emphasize an important as- 
pect of the presented discussion: IC scattering of the 
starlight emission by the synchrotron-emitting elec- 
trons is inevitable, and involves neither the unknown 
target photon field, nor the additional unconstrained 
source of the ultrarelativistic high-energy electrons. In 
particular, we 'reconstruct' the electron energy distri- 
bution from the known broad-band synchrotron spec- 
trum of a given jet region, and then estimate the IC 
flux for the known target photon field. Therefore, our 
results, presented in detail in [18j are independent on 
any model of particle acceleration. 



2. THE MODEL 

2.1. Electron energy distribution 

Detailed broad-band (radio-to- X-ray) observations 
of the M 87 large-scale jet circumscribe well the syn- 
chrotron spectrum of knot A, the brightest part of the 
jet in the radio and optical regimes. Hereafter we re- 
fer to the analysis by (2jj, which indicates that the 
energy distribution of the synchrotron-emitting elec- 
trons (as measured in the jet commoving frame de- 
noted by primes) can be approximated by a broken 
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7 p for 7 < 7b r 
7 « r7 -(f+«) for 7 > 7br 



(1) 



with the spectral indices p = 2.3 and q = 1.6. 
The electron break Lorentz factor 7b r ps 2.7 x 
10 6 <5~ - 5 BZ4' 5 corresponds to the observed syn- 
chrotron break frequency ^t> r ps 10 15 Hz for the 
emitting plasma magnetic field B = B-4 10~ 4 G and 
Doppler factor 8. The cut-off energies 7 m i n and 7m ax 
are basically unconstrained, but the synchrotron ori- 
gin of the X-ray jet emission indicates 7 roax /7br > 
10 — 100. We normalize the number of radiating elec- 
trons to the observed (isotropic) synchrotron luminos- 
ity of knot A at the break frequency ps 3 x 10 41 
erg s" 1 H3- 

2.2. Starlight photon field 

For the target starlight photons at the position of 
knot A we assume roughly isotropic distribution in 
the galactic rest frame and strongly anisotropic in 
the jet rest frame, due to the relativistic jet velocity. 
We take the characteristic observed frequency of the 
optical-NIR bump due to the elliptical host of M 87 
ps 10 14 Hz In). We also evaluate the appro- 
priate starlight energy density at the position of knot 
A directly from the observations of M 87 host galaxy, 
assuming King profile for the stellar emissivitv_ at the 
distances > 1" — 2" from the M 87 nucleus 
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for r < 7"t 
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where r is the radius as measured from the galactic 
center, r c is the core radius for the galaxy and r t is 
the appropriate tidal radius. We normalize this dis- 
tribution to the observed luminosity density profile in 
the / band 0, with the parameters r c = 0.55 kpc 
and rt = 68 kpc. We evaluate further the /-band 
stellar energy density at the position of knot A as 
C4tar(l kpc) ps 10~ 10 erg cm' 3 , which is in fact a very 
safe lower limit (see a discussion in [Ts^l. 



2.3. IC emission and kinematic factors 

We evaluate the high-energy emission of knot A 
due to IC scattering on monoenergetic and mono- 
directional (in the jet rest frame) starlight photon 
field, including properly relativistic effects in the 
Klein-Nishina regime, using the approximate expres- 
sion given by pj. We also consider the bulk Lorentz 
factor of the M 87 jet at the position of knot A, T, to 
be ~ 3-5, and the jet viewing angle 9, ~ 20° -30° @. 
For such a choice of T and 9 one gets the jet Doppler 
factor of the discussed region 8 — [T(l — (3 cos 6*)] _1 ~ 
1-3. 



3. THE RESULTS 

Spectral energy distribution of the high-energy 7- 
ray IC emission of knot A is presented in figure 1 for 
different magnetic field intensities, jet viewing angles 
9 = 30°, 20°, and the bulk Lorentz factors T = 3 
and 5. The Thomson part of this emission extends 
up to photon energies of about 10 10 — 10 11 eV. Below 
we compare the expected IC fluxes for different pho- 
ton energies to the upper limits imposed by the ob- 
servations of EGRET observatory and ground-based 
Cherenkov Telescopes: HESS, HEGRA and Whip- 
ple, and next we compare the implied lower limits on 
the jet (knot A) magnetic field with the appropriate 
equipartition value, 



B oq = 330<T 5/7 fJ,G 



(3) 



as given by Q. Note, that the adopted value of B eci 
refers to the energy equipartition between the jet mag- 
netic field and ultrarelativistic electrons, with possi- 
ble contribution from the non-radiating particles ne- 
glected. 

3.1. EGRET observations 

EGRET observations of the Virgo cluster give the 
photon flux F(> 100 MeV) < 2.2 x 10~ 8 cm- 2 
s _1 which implies B > 30 — 100 fiG correspond- 
ing roughly to B/B cq > 0.1 — 0.5, for any choice 
of the kinematic factors considered here. This con- 
straint does not necessarily imply a departure from 
the magnetic field-radiating particles energy equipar- 
tition but, interestingly enough, already excludes a 
class of models involving a very weak jet magnetic 
field. 

3.2. Whipple observations 

Whipple observations give the 99% CL upper limit 
to the VHE 7-ray photon flux of M 87 radio galaxy 
F(> 0.4 TeV) < 6.9 x 10" 12 cm" 2 s" 1 [Hfl. Indepen- 
dently on the EGRET constraints, the Whipple ob- 
servations imply thus magnetic field intensity within 
the discussed jet region B > 50 — 60 fiG, or in other 
words B/B cq > 0.2 — 0.4 again for any choice of the 
kinematic factors considered here. 

3.3. HEGRA and HESS observations 

HEGRA observations resulted in a 4<r detection of 
high-energy 7-ray flux from the direction of M 87 with 
the photon flux F(> 0.73 TeV) ps 0.96 x 10~ 12 cm" 2 
s _1 ;2]. Recent HESS observations resulted also in 
the detection of the M 87 system at the 3 — 4er level, 
however with the photon fluxes F(> 0.73 TeV) ps 0.4x 
10~ 12 cm" 2 s- 1 in 2003, and F(> 0.73 TeV) ps 0.15 x 
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Figure 1: Spectral energy distribution of high-energy 7-ray inverse-Compton emission of knot A for the jet magnetic 
field B = 0.3 x 10" 4 G, 10" 4 G, 3 x 10" 4 G, B = 10" 3 G (as indicated at each panel), the jet viewing angles 6 = 30° 
and 20° (blue and red curves, respectively) and the jet bulk Lorenz factors T = 5 and 3 (solid and dashed lines, 
respectively) . 



10~ 12 cm~ 2 s _1 in 2004 y|. This clearly indicates a 
variability of the high-energy 7-ray emission of this 
source, and therefore gives the upper limits for the 
VHE radiation of knot A. 

Figure 2 shows the expected photon flux of knot A 
at the observed photon energies hv® > 0.73 TeV as 
a function of the magnetic field intensity, for the jet 
viewing ang les 9 = 30° and 20°, and the bulk Lorentz 
factors r = 3 and 5. The vertical lines indicate the 
appropriate equipartition magnetic field. The most 
recent HESS observations imply therefore the mag- 
netic field within knot A as strong as B > 300 //G, 
i.e. B/B cq > 1 - 2. 

4. CONCLUSIONS 

Our study indicates that the magnetic field within 
the brightest knot A of M 87 jet is most likely stronger 
than the equipartition value that refers to the radiat- 
ing electrons alone, B > 300 /zG > B eq (if the jet 
viewing angle is in the range 9 = 20° — 30° and the jet 
bulk Lorentz factors T = 3 — 5). On the other hand 
the upper limit of the magnetic field intensity within 
knot A can be found from the upper limit imposed on 



the magnetic field energy flux, Lb = ^R 2 cT 2 B 2 < Lj, 
where Lj ~ few x 10 44 erg s" 1 is the total power of 
M 87 jet 01 an d R ~ 60 pc is the radius of radio 
structure at the position of knot A. This gives roughly 
B max < 1000 //G. The comoving energy density of 
the magnetic field within knot A is therefore limited 
roughly by 4 x 10~ 8 erg cm~ 3 > £/ B > 4 x 10 -9 erg 
cm -3 . 

VLBI measurements often allow one to infer mag- 
netic field intensity from the low-frequency spectral 
break in the radio emission of the (sub) pc-scale jet 
modeled in terms of synchrotron self-absorption pro- 
cess. In the case of M 87 jet, this method gives 
Svlbi < 0.2 G at Evlbi < 0.06 pc 0. If the mag- 
netic energy flux in a jet is constant, then one should 
expect magnetic field intensity at the position of knot 
A to be roughly B = (Tvlbi/T) (i? V LBi/i?) B VLBI <C 
300 /iG, where we put TyLBi/r « 2 and the jet radius 
at the position of the considered knot R = 62 pc 0] . 
Hence one can conclude that some kind of amplifica- 
tion of the jet magnetic field has to take place between 
the parsec and kiloparsec scales, although all the esti- 
mates presented above should be taken with caution, 
as some arbitrary assumptions on the jet magnetic 
field structure were invoked. 
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Figure 2: Constraints on the jet magnetic field within knot A imposed by the HEGRA and HESS observations (dotted 
horizontal lines), for the jet viewing angle 8 = 30° and 20° (blue and red curves, respectively), and the jet bulk Lorenz 
factors r = 5 and 3 (solid and dashed lines, respectively). Vertical lines denote the equipartition magnetic field for 
T = 5 and 3 (solid and dashed lines, respectively) and 8 = 30° and 20° (blue and red curves, respectively). 



We speculate that the suggested amplification of 
the jet magnetic field is due to an action of the tur- 
bulent dynamo processes discussed in this context by, 
e.g., jfj: the Kelvin-Helmholtz instabilities occurring 
inevitably at the edges of the jet are supposed to cre- 
ate large-scale eddies at the flow boundaries, which 
amplify the jet magnetic field and develop highly tur- 
bulent mixing layer between the jet and the surround- 
ing medium. Such turbulent shear layers play a crucial 
role in mass entrainment (and so deceleration) of the 
FR I outflows 4], and in acceleration of the jet par- 
ticles [laL influencing also polarization properties of 
the jetsQ. 

In a framework of this model, on the small scales 
the M 87 jet has to be dynamically dominated by 
the (cold) particles, in order to allow for the tur- 
bulent dynamo process to proceed at all. On the 
larger scales the mass entrainment process decelerates 
the flow gradually, slowly amplifying the jet magnetic 
field (to the value exceeding at some point the energy 
equipartition with the radiating electrons), and creat- 
ing a turbulent boundary layer that spreads into the 
jet interior. This process accomplishes at ~ 1 — 2 kpc 



from the active center (knot A and beyond), where the 
jet magnetic field reaches an approximate equiparti- 
tion with the energy density of the jet particles, while 
the jet flow itself becomes fully turbulent, disappear- 
ing into the outer amorphous radio lobe. 
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